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Abstract 



We investigate the correspondence between classical and quantum mechanics 
for periodically time dependent Hamiltonian systems, using the example of a 
periodically forced particle in a one-dimensional triangular well potential. In 
particular, we consider quantum mechanical Floquet states associated with 
resonances in the classical phase space. When the classical motion exhibits 
su&harmonic resonances, the corresponding Floquet states maintain the driv- 
ing field's periodicity through dynamical tunneling. This principle applies 
both to Floquet states associated with classical invariant vortex tubes sur- 
rounding stable, elliptic periodic orbits and to Floquet states that are associ- 
ated with unstable, hyperbolic periodic orbits. The triangular well model also 
poses a yet unsolved mathematical problem, related to perturbation theory for 
systems with a dense pure point spectrum. The present approximate analyt- 
ical and numerical results indicate that quantum tunneling between different 
resonance zones is of crucial importance for the question whether the driven 
triangular well has a dense point or an absolutely continuous quasienergy 
spectrum, or whether there is a transition from the one to the other. 
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I. INTRODUCTION 



Ever since Bohr's formulation of the correspondence principle ffl, the relation between 
classical and quantum mechanics has been a subject of great conceptual interest. The 
simplest type of time- independent, bounded classical Hamiltonian system is integrable, i.e., 
if there are as many (n) independent constants of motion as degrees of freedom. The phase 
space of such a system is completely stratified into invariant n-tori and the quantum 
mechanical energy eigenfunctions and eigenvalues can be determined approximately with the 
help of the Einstein-Brillouin-Keller (EBK) quantization rules If, on the other hand, 

the classical system is chaotic, then semiclassical quantization relies on classical periodic 
orbits f^-|TTl. "Typical" classical Hamiltonian systems, however, are neither completely 
regular nor completely chaotic, but exhibit an intricate coexistence of regular and apparently 
stochastic motion [TJ]. An important example of this type is the three-body Coulomb 



system. The application of both torus quantization and a periodic orbit cycle expansion 
to this system has resulted in a semiclassical understanding of the spectrum of the Helium 
atom 0. 

Much insight into the classical-quantum correspondence has been gained by the study 
of simple model systems. A popular model, the periodically kicked rotor, has led to the 
discovery of quantum mechanical suppression of classical diffusion [14 , a phenomenon that 



could subsequently be explained by a mechanism closely related to Anderson localization in 
disordered one- dimensional crystal lattices []l5i|l6| . For a general quantum system defined 



by a periodically time dependent Hamiltonian operator H(t) = H(t + T) the role of the 
stationary states is taken over by the Floquet states Hl7|p!8| . If the operator 



H := H{t) - ihd t , (1) 



acting on an extended Hilbert space of T-periodic functions |19| , has T-periodic eigenfunc- 
tions, i.e., if the eigenvalue equation 

\H{t) - ihd t ]u a (t) (t) (2) 



can be solved with periodic boundary conditions in time, u a (t) = u a (t+T), then the Floquet 

states 

i>a(t) '■= u a (t) exp(— ie a t/h) (3) 

are solutions to the time-dependent Schrodinger equation. The eigenvalue equation (0) now 
plays a role analogous to that of the stationary Schrodinger equation for time-independent 
systems, and the objective of a semiclassical theory is the approximate computation of the 
Floquet eigenfunctions u a (t) and the quasienergies e a , starting again from invariant objects 
in the classical phase space. 

However, periodically time dependent quantum systems pose some mathematical diffi- 
culties that are rarely met in the time-independent case. Consider a Hamiltonian operator 
of the form H(t) = H + (5Hi(t), where only H\(t) = Hx(t + T) carries the time dependence, 
and /3 is a dimensionless coupling constant. Assume further that H possesses a discrete 
spectrum of energy eigenvalues E n (n = 1,2,..., oo), with corresponding eigenfunctions (p n . 
Then for j3 = the Floquet states can be written as 

^(n,m)(t) = (W m "') exp[-i(E n + mhw)t/h] (4) 

with oo = 2tt/T. If m is an integer number, the Floquet function U(n tm \(t) = ip n e muJt is T- 
periodic, as required; the corresponding quasienergies are £( n ,m) = E n + mhu. The index a 
in (0) thus becomes a double index, a = (n,m), and the quasienergy spectrum is given by 
the energy eigenvalues modulo %u. Therefore it can be divided into Brillouin zones of width 

hoj. 

Apart from special cases, this spectrum will be dense on the real axis. But what happens 
when (3 > 0? Rigorous perturbation theory for such systems with a dense point spectrum is 
anything but trivial [^,^1|] . The conditions for the perturbed system to retain a pure point 
spectrum constitute a crucial unsolved problem. Partial answers to this problem of "quan- 
tum stability" |TL| have been obtained only recently: it has been shown by Howland |22 



that even for arbitrarily large /3 the quasienergy spectrum remains a dense pure point spec- 
trum, provided the gaps E n+ i — E n between the energy eigenvalues of H grow sufficiently 



rapidly with n, and the perturbation H\(t) is bounded. This result applies, for instance, to a 
periodically forced particle in a one-dimensional box. The condition on the boundedness of 
the perturbation can be relaxed in the case of certain other periodically forced anharmonic 
oscillators |[23|| , but the condition on the growth of the energy gaps between the eigenvalues 
of H remained essential also in subsequent studies • 

A model that violates both propositions of Howland's theorem is given by a periodically 
driven particle in a triangular well potential: the unperturbed system has a pure point 
spectrum, with the gaps between successive energy eigenvalues E n decreasing asymptotically 
as n -1 / 3 , and its dipole operator is unbounded. This model emerges, for example, in the 
description of far-infrared irradiated semiconductor hetero junctions |26|, or of microwave- 



driven charge carriers in superffuid helium [27|. 



This model is the subject of the present work. Although fairly simple, it provides in- 
teresting insights into the classical-quantum correspondence for periodically time- dependent 
systems with a mixed classical phase space. As we will demonstrate, many features of the 
quantum system have a very close link to its classical counterpart, both in the classically 
regular (elliptic) and chaotic (hyperbolic) case. Nevertheless, quantum mechanical tunneling 
is of crucial importance to understand prominent features of both the quasienergy spectrum 
and the Floquet eigenfunctions. 

The tunnel effect might also play a central role in the question whether there is a pure 
point spectrum or an absolutely continuous one. From a rigorous point of view this difficult 
question is still open, but it has been argued that the quasienergy spectrum of the peri- 
odically driven triangular well exhibits a transition from a point spectrum to a continuous 



spectrum at a nonzero, finite value of the driving amplitude p8| , p9| . Although there are 
well studied examples of perturbations through which a discrete spectrum becomes contin- 
uous 



50[j , a spectral transition in a system like the driven triangular well would be quite 
unusual, and therefore interesting. If there is such a transition, there must be a mechanism 
producing it; however, what could that be? 

Thus, the aim of this work is twofold. First, we wish to provide instructive examples 



of classical-quantum correspondence and non- correspondence. Due to its simplicity, the 
triangular well model might even be useful for courses in elementary quantum mechanics, 
despite the open mathematical problems. Second, we wish to stimulate further rigorous 
research on the nature of the quasienergy spectrum of periodically time dependent quantum 
systems. Previous work P0|-f25|] has addressed this problem in a quite formal way. However, 
a physically oriented starting point, such as a semiclassical analysis, might also prove fruitful, 
and provide complementary information. 

We proceed as follows. Section 2 introduces the classical model and states the results of 
a standard resonance analysis. Unlike the periodically forced particle in a box considered by 
Reichl and Lin pl| , |32| , the triangular well model has the property that even for arbitrarily 
small driving amplitude there exist arbitrarily large primary resonances in the classical phase 
space. Therefore, even for small amplitudes an infinite number of quantum mechanical 
Floquet states will be influenced by these resonances. Section 3 contains a semiclassically- 
motivated description of near-resonant Floquet states, followed in Section 4 by a more refined 



analysis, which is an extension of previous work by Berman and Zaslavsky ]3"3" |. Section 5 
compares the analytically calculated quasienergy spectrum with numerical data. The central 
Section 6 explores the influence of classical phase space structures on the Floquet states. 
We do not focus on questions pertaining to the resonance overlap |pli|32| , ^4| , |35| , but rather 
try to stretch the single-resonance picture as far as possible. An interesting result that 
will be obtained in this way is the existence of a closed link between states associated with 
regular classical dynamics and others linked with apparently chaotic parts of the classical 
phase space: Floquet states associated with unstable, hyperbolic periodic orbits appear as 
certain excitations of states most closely associated with invariant manifolds surrounding 
stable elliptic periodic orbits. The final discussion then ties together some loose ends, and 
returns to the question whether or not there could be a spectral transition. 
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II. THE CLASSICAL TRIANGULAR WELL 



A classical particle of mass m that is confined to a one-dimensional triangular well 
and driven by a monochromatic force of amplitude F and frequency us is described by the 
Hamiltonian 

P 2 

H(P,X,t) = — + V a (X) + FXcos(cot) (5) 

with V a (X) = aX for X > and V a (X) = oo for X < 0. Introducing dimensionless 
variables, 

x = muj 2 X/a 
p = ujPja 

T = Ult (6) 

this Hamiltonian can be written in the form H(P,X,t) = a 2 /(muj 2 )H(p,x,T), where the 
scaled Hamiltonian 

p 2 

H(p,x,r) — — + Vi(x) + /?xcos(r) = H (p,x) + /5xcos(r) (7) 

depends only on the scaled amplitude (3 = F/a. A detailed analysis of the "kicked" version 
of this system, where the sinusoidal time dependence is replaced by a sequence of 5 kicks 
with alternating signs, has been presented by Shimshoni and Smilansky |27 • 



To express H(p,x, r) in terms of the action- angle variables of the undriven system defined 
by H (p,x), we first exploit the relation between the action / and the energy E: 

I = - [ E ^2(E-x)dx = ^{2Ef 12 ■ (8) 

TT JO 07T 

Hence, 

H (I) = (3nI) 2 / 3 /2 . (9) 

The oscillation frequency of the undriven particle in the triangular well is then found 
from Hamilton's equation 



and the angle variable is 9 = Q(I)t, restricted to the interval from to 2n. 

Imposing the initial condition x(r = 0) = 0, the unperturbed trajectories can be written 

as 

W 3/2/3 L 02 ^ 
x{0) 



(Svr) 1 / 3 V n =l ™ 
Thus, the full Hamiltonian becomes 



J2/3 / oo fi \ 

^-E^cos^) . (11) 



/o r\2/3 o /q r\2/3 oo -i 

jf (I, 0, r) = 1 J + ^(3vr/) 2 / 3 cos(r) - (3 { [ £ - [cos(r + n0) + cos(r - nB)\ . 

(12) 

Primary resonances occur when the particle's action is such that the cycle time of unper- 



turbed motion coincides with an integer multiple, N, of the period of the driving force ||12| , |36 
Then the argument r — n9 in the above expansion becomes stationary for n = N . The res- 
onant action I N can thus be found from the equation $1(1 n) = l/N. By fllUp, 

n 2 N 3 , . 

Is = — ■ (13) 

Approximating Hq(I) quadratically around In and keeping only the resonant term of 



the perturbation (taken at I — In), one obtains the familiar pendulum Hamiltonian [|T2|j36 
H(I, 9, r) « H (I N ) + n(I N )(I -Is) + 1^(1- In) 2 ~ (3 cos(r - N6) , (14) 

where 

1 



M _1 ^ dm 



dl 



T n 2 N 4 



(15) 



This Hamiltonian describes the regular, resonant motion close to In- Within the pendulum 
approximation, the half-width of the resonances is 
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AI N = 2irNyP . (16) 
According to the Chirikov criterion p~2L|36f, the N-th resonance overlaps with the (N + l)-th 



when Jjv+i — In = A/at+i + Aijv- The corresponding critical values of the driving strength 
are 

Pn ~ 16 \N* + N + 1/2 ) ' 1 ] 

For example, /?( c) = 0.537 and (3^ = 0.586. 

Summing up, the model of the classical driven particle in a triangular well has the 
following three essential features: (i) there are an infinite number of primary resonances, 
with the iV-th resonant action 1^ being proportional to N 3 ; (ii) the half- widths Aijv grow 
quadratically with the order N of resonance; (iii) within the pendulum approximation, all 
resonance overlaps occur in the interval 0.537 < (3 < 0.617 of the control parameter f3. 

Of course, the pendulum approximation can only give a coarse description of the actual 
dynamics, and is restricted to moderate values of (3. The larger (3, the more important 
are the stochastic layers surrounding the resonances JT2J. When (3 is increased beyond the 



perturbative regime, the widths of the zones of resonant, mainly regular motion in phase 
space do no longer grow proportional to ^/J3, as suggested by fllfip, but they actually shrink. 
The \f/3 growth is overcompensated by the growth of the stochastic layers. This fact is 
illustrated by the following two figures. Fig. 1 shows a Poincare surface of section for the 
driven particle in a triangular well with f3 = 0.5, taken for r = 37r/2 mod 2ir. Apart from 
small secondary resonances, the most prominent features are the island of regular motion 
for iV = 1 in the lower left corner, and the two islands for N = 2. The stochastic layers, 
which for smaller (3 merely cover limited areas bordering the individual resonance zones, have 
already merged into a connected sea of stochastic motion. The locations of the resonant 
islands agree reasonably well with the predictions of (|T3"D, but for N = 2 the half-width is 
clearly smaller than A/2, as given by ([II]). As shown in Fig. 2, for (3 = 1.0 the islands are 
squeezed to a small fraction of their previous size. But still, there exists regular motion 
"inside" the resonant islands, although (3 is significantly larger than (3^ . 
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III. SEMICLASSICAL APPROACH 



The scaled classical Hamiltonian (0) depends only on a single dimensionless parameter. 
A second one appears in quantum mechanics: Planck's constant Ti is scaled to 

h eff := mu 3 h/a 2 (18) 

as a consequence of the transformations (||). The model of the periodically driven triangular 
well, therefore, is of paradigmatic simplicity for a study of the classical-quantum correspon- 
dence: one parameter, (3, controlls the degree of nonintegrability in the classical system, 
the other, h e ff, controlls the scale at which its quantum mechanical counterpart can resolve 
classical phase space structures. 

Since I — In and 6 constitute a canonically conjugate pair, the approximate pendulum 
Hamiltonian ( |14D can be quantized by the replacement 

< 19 > 

This is essentially a semiclassical argument, since it ignores the fact that the canonical 
transformation (p, x) — > (I, 6) does not commute with the quantization operation. The 
resulting Schrodinger equation 



d 



H (I N ) + Wn)^^ - -[3cos(N6 - r) 



i 36 2M 09 2 
can be solved with the help of the ansatz 



W,t) (20) 



^(0, r) = X (z) e W (-iWr/fi eff ) , (21) 

where the new variable z is given by 

z = (N6-t)/2. (22) 

Inserting (|21|) into (|20| ) and using the resonance condition Q(In) = one obtains the 

Mathieu equation 



9 



'A 2 

dz 



2 + a-2qcos(2z) xO) = 



(23) 



with parameters 



A\M\ 



0- 



(24) 
(25) 



The general solution to this equation has the Floquet form x( z ) = Pv{z)e tuz , with a charac- 
teristic exponent v and a 7r-periodic function P v {z) = P v (z + tt). Thus, the corresponding 
wave functions ip(9,r) have the form 



4>(0,t) 



(N9-t\ ( NO-t Wt 
Pv — exp iv i- 



h 



eff. 



(26) 



Since 9 is an angle variable, ip must be 27r-periodic in 9. This requirement now restricts the 
possible values of v to the discrete set 

2j 



v = 



N 



j = 0,l,2,...,iV-l. 



(27) 



According to (p5|), the parameter q entering the Mathieu equation is fixed by the choice of 
the amplitude (3. The allowed values u(j) can occur as characteristic exponents of solutions 
to the Mathieu equation only for certain values (labelled by the integer k) ak{v{j)-, q) of the 
other parameter a. These values q) finally determine, by (E¥), the allowed Wj t k- 

The wave functions (EH) are now fully specified. Since they already are of the required form 
"0j,fc(^) r ) = u j,k(9,r) exp(— ie^fcT /h e ff), with Uj t k(9,r) = Uj t k(9, t + 2n), the quasienergies 
for the approximate equation (|2~0D have been found: 



h 



j. A; 



mod h 



eff 



k 2 ff N 2 7 ' 



mod h 



eff 



The cases N = 1 and N = 2 are of particular interest. For N = 1 there is j 



(28) 







only, and the a^(0, q) coincide with the well known characteristic values |37| that give rise 



10 



to 7r-periodic Mathieu functions. Those values associated with even Mathieu functions are 
customarily denoted |57[] as a (q), a 2 (g), a 4 (g), . . ., and those associated with odd functions 
as 62(5), 64(9)) be(q), ■ ■ ■ With the definition 



«o,fc(g) := 



Ofc(g) , k = 0,2,4,... 
h+i(q) , fc = 1,3,5, . . . 

the approximate quasienergies for an N = 1 resonance can be written as 



(29) 



H (I. 



h 



N , 



eff 



\\M\ 



a ,k(q) mod h ef f 



(30) 



For N = 2 there are two groups of states, labelled by j = and j = 1. Those with j = 1 
have to be constructed from the 27r-periodic Mathieu functions. The required characteristic 
values ai t k(q) = a k (l,q) are 



b k+l {q) , k = 0,2,4, 



a k {q) 



(31) 



k — 1 . 3. 5 , 



Thus, the quasienergies are given by 



Hod 



mod /i e /j 



(32) 



for j = 0, 1 and fc = 0, 1, 2, . . . 

This spectrum has a simple interpretation 



The Mathieu equation ( p3|) can be 



regarded as a stationary Schrodinger equation for a fictitious particle moving in a cosine 
lattice, and the solution for the iV-th resonance requires periodic boundary conditions after 
iV cosine wells. Thus, N = 1 leads to a single-well potential, and iV = 2 to a double 
well. The quasi-energies of the original problem correspond to the energies of the fictitious 
particle. The eigenstates in the single well are simply labelled by the quantum number 
k. The spectrum of the double well, however, consists of a sequence of doublets, each one 
associated with a symmetric and an antisymmetric eigenfunction. The absolute value of the 
difference £1^ — eo,k — heff/2 =: (Ae) k of eigenvalues within each doublet corresponds to the 
tunnel splitting. The well known asymptotic behavior of the characteristic values [j37j]. 
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[77 k/2+3/4 

b k+ M - a k (q) ~ 2 4fc+5 y-^— exp(-4 v ^) , (33) 

shows that the tunnel splitting becomes exponentially small with ^fq, i.e., with the square 
root of the driving amplitude (3 (see (p5D). That is easy to understand: g determines the 
depth of the individual cosine wells; the larger q, the weaker the tunneling between them. 
In particular, the ground state splitting, i.e., the (quasi-) energy difference between the two 
states most tightly bound by the double well potential, is 

(Ae)o = |^ (h(q) - oo(ff)) 

yV 4 exp(-V^) (34) 

with g = IGir 2 (3 /hlj f . It is quite remarkable that the approximate construction of near- 
resonant Floquet states for N > 1 accounts for the tunnel effect, even within the present 
semiclassical approach. 



IV. QUANTUM MECHANICAL APPROACH 

The considerations in the previous section, albeit instructive, have several flaws. Most 
of all, it is not easy to control the accuracy of the semiclassical replacement (|l9l) , and an 
expression like H (In) will, in general, not correspond directly to a quantum mechanical 
energy eigenvalue. It is, therefore, desirable to analyze the role of classical resonances in 
quantum systems in a strictly quantum mechanical way, without direct recourse to classical 
mechanics. In this section, we sketch such a quantum mechanical approach. It is basically a 
generalization of the previous work by Berman and Zaslavsky [[£J . Although it also confirms 
the main conclusions of the semiclassical computation, its primary merit lies in the fact that 
it clearly shows the physical and mathematical difficulties that will be encountered in the 
attempt to improve the approximation. 

We consider the Hamiltonian fl7|) and denote the unperturbed energy eigenstates and 
eigenvalues of the triangular well by \n) and E n , respectively. We assume that the r-th 

12 



eigenstate is "resonant", i.e., that the level spacing close to the r-th eigenstate is approxi- 
mately equal to h e ff/N for some integer N: 



, _ dE n 



dn 



heff 



(35) 



N 

r 

This is the quantum mechanical equivalent of the classical resonance condition £1(1 n) = 1/N. 
Only for special choices of h e ff will it be possible to obtain an exact identity E' r = K e ff/N\ 
in the general case we choose r such that the absolute value of the difference 

8:=E>- % -f (36) 

is minimized. At most, 5 will be of the order of E". We also stipulate that r be so large 
that for levels close to E r the spacing between E n and E n+ ^ varies only slowly on the scale 
of h e ff, which implies the inequality 

K < K ■ ( 37 ) 

We are interested in a solution \ip(r)) to the time-dependent Schrodinger equation that 
consists mainly of a superposition of near-resonant eigenstates. Thus, a reasonable ansatz 
is provided by 



(r))=E c n(r)\n)e^Li(E r + (n-r)^j^} , (38) 



eff 



which immediately leads to 

n — r 



ih e f f c n (r) = [E n - E r - h e ff—— j c n (r) + /5cos(r) ^(n|x|m)c m (r) exp[z(n -m)r/N] . 

(39) 



Keeping in mind the inequality (|37|) , we expand the eigenvalues quadratically around E r , 

E n ^ E r + {n-r)E' r + hn-r) 2 E'; . (40) 

According to our propositions, the time-dependence of the coefficients c n (r) should be 
weak, compared to that of the exponentials exp(i£r/N) with a nonzero integer £. Corre- 
spondingly, we keep only the stationary terms m = n ± N on the r.h.s. of (0). Thus each 
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coefficient c n (r) is coupled only to c h ±n{t), c 71 ±2Jv(t), and so forth. One obtains AT disjoint 
sets of coefficients, which we denote by Within each set, j — 0, 1, . . . , TV — 1 is constant, 
and n is restricted to the values n = r + j + mN, where m is an integer. In addition, we 
assume that the dipole matrix elements (5{n\x\n ± N) can be approximated by a constant 
V, independent of n. Using now the new index m to label the coefficients, one arrives at 
the N decoupled sets of equations: 



with Wj still to be determined. Assuming that fj(<p) is 27r-periodic, the auxiliary function 
gj(ip) = fj(<f)exp(ij<f/N) is 2A r 7r-periodic. The resulting equation for gj((p), i.e., 



iteff&Hr) = \N 2 E';{m + j/Nfc^{r) + N5(m + j /N)c^(r) + V - (c^ +l {r) + d 




(41) 



To solve these equations, we set 




(42) 



( 



1 9 „ d 2 NS d , « \ , N 

-2 NE "w + -T v + Vms ^- w i)*^ = ' 



(43) 



can then be transformed into a standard Mathieu equation; the ansatz 



g^) = eM-^zS/(NE';)] X ^) 



(44) 



with ip = 2z leads to 




(45) 



with the parameters 




(46) 



AV 



(47) 



Q = 



N 2 E'J ' 
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Starting from a general Floquet solution to the Mathieu equation (f45|), Xj( z ) — 
P^(j)(z) exp(iu(j)z) , the function gj((p) can now be expressed as 

9j (<p) = exp [i^ - iJEL^J . (48) 

Since, on the other hand, gj(f) = fj(^>) exp(ij(p/N), one finds fj(<p) = P v (j)(ip/2), together 
with 

u(j) = ^ + ; j = 1, . . . , N - 1 . (49) 

As in the previous section, these characteristic exponents admit only a discrete set of allowed 
values for the parameter a; these values are again denoted by a k (i/(j), q). By (|46|), they lead 
to a discrete set of possible values Wj k for the quantities Wj. Now the expansion coefficients 
are completely specified: 

1 r 2n 



cliHr) 



2tt Jo 

(m) 



/ d<p P v (j),k(<p/2) exp(-zm^ - iW jt k T/h eff ) 

J 



= xTk >e M-iW jik T/h eff ) . (50) 

Thus, the resulting approximate Floquet solutions \i>j,k{ T )) to the Schrodinger equation can 
finally be written as 

foOO) = E 4?k + J + mN)e~^ exp j-i fs r + W hk + ft, - — 



l eff 



NJ h ef fl 

Uj, k (r)) exp(-ie j>k T/h ef f) . (51) 



Using (|46|), the quasienergies e^* = _E r + + h e ffj/N (mod ft e //) are 

mod H e ff . (52) 
The structure of the near-resonant quasienergy spectrum obtained in this way is equal to 



£j,k 



the structure of (28), but the parameters differ. These differences are small in the regime 
where the approximations are valid. First, the exact energy eigenvalues for the triangular 
well potential are given by 

(h 2 \ 1/3 

E »=[^j L ) z m (53) 
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where the numbers z n (n = 1,2,3,...) denote the negative zeros of the Airy function 



Ai(z) |7]: 



(54) 



with 



/(*W /3 (i + ^- 2 -^ 4 + -- 



5 

36 



(55) 



If this expansion is approximated by its leading term, and if the Maslov correction resulting 
from the right turning point of the motion in the triangular well potential is neglected, one 
has 



and hence 



y(--^)) 2/3 -(3W2) 2/3 



E n ^-(3nh eff nf 3 



(56) 



(57) 



Since the quantum number r of the resonant state is large by assumption, these "semiclas- 
sical" approximations are appropriate. Thus, within these approximations the eigenvalues 
E n agree with those obtained from the classical Hamiltonian Hq(I) by the Bohr-Sommerfeld 
replacement I — > h e ff n , see (f|). With I N « ^ e // r , the parameter M _1 in (|28|) approximates 
E';/hl ff in flU), see P). 

The dipole matrix elements of the triangular well can be computed analytically fl39| : 



(n\x\m) 



(2h eff ) 2 / 3 z n /3 



n = m 



-(2h eff ) 2/3 /(z n - z m f , n/m 



(58) 



Using the arithmetic mean 







V — — {{r\x\r + N) + (r[a;|r — N)) 



(59) 



for the coupling constant V appearing in (f4"Tp, and expanding z r ±N to first order in N/r, 
one finds 
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V&-/3. (60) 
Then the semiclassical and the quantum expression for the Mathieu parameter q, ( |25|) 



and fl47D, coincide. 

The characteristic exponents (|27| ) are exactly equal to those in P9| ) only if 5 vanishes, i.e., 
if E' r is exactly equal to h e ff/N. In the general case, there can therefore be a slight difference 
between the values of the parameters afc(z/(j),g) in both approximation schemes. However, 
since these parameters lie between the characteristic values ak(q) and bk + i(q), and since 
bk+x(q) ~ a k{q) vanishes for large q according to (|55D, this difference is visible only for small 
q. The term —6 2 /(2E f 1 !) appearing in (|52j) makes sure that the approximate quasienergies 
merge into the quadratically approximated eigenvalues ( f40"l) for vanishing q. Since 5 is, at 
most, of the order of E", one sees from (|37| ) and ( j36|) that this term is small compared to 
heff/N. Thus, there is a one-to-one correspondence between the semiclassically calculated 
spectrum (^) and (|52|). 

The approximations made in the quantum calculation have their obvious analogues in 
the previous section. There is one detail that deserves particular attention: it has tacitly 
been assumed that the Fourier coefficients Xjj? °f the approximate Floquet states (|5T| ) decay 
so rapidly that the non-existing "eigenstates below the ground state carry zero weight. 
The approximations leading to the Mathieu equation introduce an additional symmetry 
(the translational invariance of the cosine lattice) which renders the problem integrable. 
The larger q, the broader the distribution of the Fourier coefficients. When the distribution 
becomes so broad that the wave function can feel deviations from the lattice symmetry, the 
Mathieu approximation becomes insufficient. 



V. THE QUASIENERGY SPECTRUM: ANALYTICAL VERSUS NUMERICAL 

RESULTS 

To compute quasienergies and Floquet states for the periodically forced triangular 
well (0) numerically, we expand the wave functions in a truncated basis of unperturbed 
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eigenf unctions (x\n) = <f n (x), with n = 1,... n max . Within this basis set, we compute 
the monodromy operator (i.e., the "one-cycle evolution operator") U(2n, 0). Its eigenvalues 
are the Floquet multipliers exp(—i27ie n /h e ff); the eigenvectors yield the Floquet functions 
u n (x, 0) at time r = 0. The numerical stability of the results obtained in this way has been 
checked by varying the truncation border n max . However, even though the results appear 
stable, it is by no means certain whether such a numerical procedure can, for sufficiently 
large basis size, actually "converge" in the strict mathematical sense. This question is related 
to a deeper physical problem which we shall discuss in the final section; for the moment, we 
take the numerical results at face value. 

The following two figures show parts of the quasienergy spectrum versus the driving 
amplitude (3. Because of the modh e ff structure of the spectrum, each of the (infinitely 
many) Floquet states has a representative quasienergy within the first quasienergy Brillouin 
zone —h e ff/2 < e < +h e ff/2. That brings about the necessity to select from all computed 
quasienergies only a small fraction for graphical display. The computed spectra contain a 
plethora of avoided crossings, and it is not always possible to select the quasienergies in such 
a way that the displayed eigenvalues appear to vary smoothly with (3. Therefore, abruptly 
ending or beginning lines visible in the figures are not caused by some sort of numerical 
instability; such discontinuities are merely unavoidable in the selection procedure. 

First, we choose h e ff = 1.33325. With this choice, one has E' 2Q = h e ff/2, so that the 
energy eigenstate r = 20 is exactly resonant for N = 2. Fig. 3 shows only quasienergies 
originating from states close to the resonant one. From (|3^) one expects two almost iden- 
tical groups of eigenvalues, separated by h e ff/2. This expectation is fully confirmed by 
the numerical result; quasienergies belonging to one of the two groups are labelled by the 
Mathieu quantum numbers k in the right margin of Fig. 3. In view of the corresponding 
classical mechanics, the fact that for (3 = 0.5 at least 2x8 states are accessible to the Math- 
ieu approximation is quite remarkable: the upper right corner of the phase space portrait 
displayed in Fig. 1 shows the area 2nh e ff, corresponding to the area "occupied" by a single 
Floquet state. Judging by the total areas of the islands of regular motion, one might expect 
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the classical N = 2 resonance to influence roughly 2x4 Floquet states, merely half the 
actual number. However, it must be borne in mind that the actual regime of influence of 
a classical resonance is not restricted to the regular islands. In their vicinity there are still 
lots of fragments of invariant manifolds, such as Cantori. 

Fig. 4 displays quasienergies for h e ff = 0.66219, so that r = 40 is exactly resonant 
for N = 2, and the number of states accessible to the Mathieu approximation is doubled. 
To check the accuracy of the previous analytical considerations, Fig. 5 shows the evalua- 
tion of (|52f) for the same parameters as used in Fig. 4, and k = 0, . . . ,9. The semiclassical 
expression M _1 here differs by merely 0.047% from E'l /h 2 e ff] the Mathieu parameter q is ap- 
proximately given by 16ir 2 (3 /ti^ff- The characteristic values q) have been computed 
by an algorithm suggested in which has been found to be both efficient and accurate. 



The agreement between Figs. 4 and 5 is quite good, although one of the approximations 
used in the derivation of fl52|) is certainly critical: the assumption that the matrix elements 
(n|x|n ± 2) close to n = r = 40 can simply be approximated by a constant is investigated 
in Fig. 6. Although these matrix elements vary substantially in the range from n = 30 to 
50, the comparison of Figs. 4 and 5 indicates that the somewhat crude approximation is 
still viable. This observation is rendered plausible by a perturbative argument: a better 
approximation to the "nearest neighbor couplings" in ([y]) would be given by a linear fit, 
j3(n\x\n ± N) ~ V + (n — r)V, leading to a perturbation 

F pert «^^(n-r)(|n)(n + iV| + |n)(n-iV|) . (61) 

The matrix elements {uj^{T)\Hp ert \uj t k{.T)) , taken with the approximate Floquet states flST]), 
then contain only contributions proportional to either e~ lT or to e +lT . Therefore, to first 
order in H pert the approximate quasienergies (|5^) remain unchanged: 

I-21T 

/ dr (u jik (T)\H pert \u jik (r)} = . (62) 
Jo 

However, the significance of the Rayleigh-Schrodinger perturbation series for operators act- 
ing in the extended Hilbert space of 27r-periodic functions |1| is questionable. This is again 

19 



a problem with perturbing of a dense point spectrum; the mod h e ff structure gives rise to 
an infinite number of small denominators in the expression for the second order perturbative 
contribution of H pert . 



VI. CLASSICAL-QUANTUM CORRESPONDENCE 

Figs. 7 and 8 show contour plots of the probability density \uj t k{x, r)\ 2 of a (numerically 
computed) member of the ground state doublet (k = 0) for the N = 2 resonance and (3 = 0.5, 
with h e ff = 1.33325 and 0.66219, respectively. In these and the following figures, the spatial 
coordinate x ranges from x min = to x max = 25, and the displayed interval of time is 8n, 
corresponding to 4 cycles of the driving force. 

The structure of these Floquet states, as well as that of the quasienergy spectrum of 
near-resonant states, can intuitively be understood if the quantum dynamics is connected 
in a direct way to the corresponding classical mechanics. Whereas quantum mechanical 
energy eigenstates of systems with a (partially) integrable classical counterpart are associated 
with invariant manifolds in the usual phase space spanned by momentum and position 
coordinates [^TJ, Floquet states of a (partially) integrable system with a Hamiltonian that 



is 27r-periodic in r are associated with invariant, 27r-periodic vortex tubes || in the odd- 
dimensional, extended phase space {(p,x,r)} Consequently, there are semiclassical 
quantization rules for Floquet states that resemble closely the well-known Einstein-Brillouin- 
Keller rules |3|-§| for energy eigenstates. In the present case of a system with merely one 
spatial degree of freedom and an effective Planck constant h e ff, these rules can be formulated 
as follows: First, the correct "quantized" vortex tubes are selected by 

j> pdx = 2-irheff (k + T^j , (63) 

where the quantization path 71 winds once around a 27r-periodic tube in a plane r = const., 
and k = 0, 1, 2, ... is an integer. The quasienergies then follow from the second rule 

£ — ——— [ {pdx — Hdr) + mh e ff , (64) 
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where the path 72 is 27r-periodic in r and lies on such a tube, with the integration extending 
over one period. The integer m accounts for the modh e ff structure of the quasienergy 
spectrum |42|j . 

It is now crucial to realize that the direct computation of Floquet states by these rules 
is possible only if there actually are vortex tubes that inherit the 27r-periodicity of the 
underlying Hamiltonian. The rules can, therefore, be applied to vortex tubes that are merely 
perturbative (and, hence, 27r-periodic) deformations of energy manifolds of the undriven 
system, or to the vortex tubes of an N = 1 resonance. Such a resonance emerges if the time 
of one cycle of unperturbed motion is equal to the driving period, 2tc. Hence the central 
elliptic periodic orbit of such a resonance is also 27r-periodic, and the rules (|63|), fl64| ) refer 
to the vortex tubes surrounding it. The fact that these tubes are not simply deformations 
of energy manifolds of the undriven system is reflected by the emergence of a new quantum 
number, k. According to (^f), the resonant ground state, k = 0, is associated with the 
innermost quantized tube surrounding the elliptic periodic orbit; the first excited state, 
k — 1, with the next largest quantized tube, and so forth. Naively, one might expect the 
hierarchy of resonant states to extend as far as there are still preserved tubes. 

Obviously, the quantum number k that appears in the semiclassical rule ( |SBD can be 
identified with the Mathieu quantum number k in (|30|) ; for example, the innermost quantized 
vortex tube corresponds to the state that is most tightly bound in the cosine potential. In 
the example shown in Fig. 1, the closed curves seen inside the N = 1 island in the lower 
left corner are sections of 27r-periodic vortex tubes with the plane r = 37r/2, and ( |63"D 
applies, provided that h e ff is small enough so that the required tubes actually "fit inside" 
the resonance. 



An N = 2 resonance, however, requires a different reasoning ||38|| . In this case, the cen- 
tral elliptic periodic orbits, as well as the surrounding vortex tubes, are 47r-periodic. Of 
course, one can still use fl63|) to select "quantized" vortex tubes and compute the associated 



semiclassical wave functions. In principle, such a procedure can be carried through in some 



detail ||43|| , but the wave functions constructed in this way are not Floquet states. Generally, 
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the association of a wave function with an invariant classical object living in the extended 
phase space is witnessed by the concentration of the probability density along the projec- 
tion of that object to the "configuration space" {(x,r)}. Hence, wave functions associated 
with a 47r-periodic vortex tube would also be 47r-periodic in time, whereas a Floquet state 
necessarily has to be 27r-periodic. 

This difficulty can only be resolved if a Floquet state is not associated with a single, 4it- 
periodic tube, but rather with both parts that result from its projection to the fundamental 
piece {(p, x, r),0 < r < 2n} of extended phase space |38fl . In this way, 27r-periodicity 



can be restored. The existence of two equivalent families of classical 47r-periodic tubes, 
displaced from each other in time by a single period, imply that the associated quantum 
mechanical Floquet states must appear in pairs. That corresponds precisely to the Mathieu 
approximation for N = 2; the eigenfunctions of the double cosine well are not localized in 
the individual wells (as the corresponding Floquet states of the original problem are not 
associated with a single invariant vortex tube), but are equally extended over both wells 
(as the Floquet states are equally associated with both tubes), so that doublets of states 
emerge, separated in (quasi-) energy by the tunnel splitting. 

The additional quasienergy separation of h e ff/2 between the two groups of states with 
j = and j = 1 in (^) has a simple explanation in terms of an analogue borrowed from 
solid state physics. If the lattice constant of a one-dimensional lattice is doubled from d to 
2d, the width of the crystal momentum Brillouin zone shrinks from 2k jd to ir/d. A similar 
"dimerization" occurs here: in the vicinity of an N = 2 resonance the classical invariant 
manifolds are 2 • 27r-periodic, so that the width of the corresponding quasienergy Brillouin 
zone should be h e ff/2, instead of fr e //- Thus, the separation by Ti e ff/2 can be regarded as 
the manifestation of a purely classical effect. In contrast, the additional tunnel splitting is 
of genuinely quantum mechanical origin. 

The "ground state" Floquet states displayed in Figs. 7 and 8 likewise both reflect the 
dynamics of the corresponding classical system and exhibit a purely quantum mechanical 
effect. A particle following the stable elliptic periodic orbit of the N = 2 resonance bounces 
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against the wall of the triangular well, is reflected, runs against the linear slope, is deceler- 
ated, reverses its direction and hits the wall again precisely Att after the previous bounce. 
There is an equivalent periodic orbit displaced in time from the first one by 2tt. The two 
quantized vortex tubes with k = surround these orbits in the extended phase space, and, 
as seen in Figs. 7 and 8, the quantum mechanical probability density of an associated Flo- 
quet state is concentrated along the projection of both tubes to the (x, r) plane. Now the 
projection of a tube that describes classical particles moving away from the wall intersects, 
at certain moments, the projection of a tube describing motion towards the wall. In the 
intersection regions there is a strong interference of both quantum mechanical "possibili- 
ties" . The two equivalent classical tubes are isolated in the extended phase space, but the 
associated wave functions communicate. 

If one forms even or odd linear combinations of the eigenfunctions of a symmetric double 
well, one obtains functions that are localized in the individual wells. In complete analogy, it is 
possible to form superpositions of the two members of the resonant ground state doublet that 
are localized along the projection of only one of its two tubes. This is demonstrated in Fig. 9, 
again for h e ff = 0.66219 and (5 = 0.5. The wave function displayed there appears to follow 
only one of the two equivalent tubes and, thus, manifestly exhibits subharmonic motion. 
However, it remains coupled to the other tube by quantum tunneling: after the tunneling 
time r tunne i := h e ff7r/(Ae) the density would be concentrated along the projection of the 
other tube, and then start to tunnel back. The approximation (|34| ) yields (Ae) ~ 8.5 -10~ 24 , 
and thus a tunneling time T tunne i ~ 4 • 10 22 • 27r. 

The same principle holds for the excited resonant states, i.e., for those states with higher 
Mathieu quantum number k. For example, Fig. 10 shows the density of a pure Floquet 
state with k = 6, for the same parameters as before. According to (|63"D, its associated tubes 
are much wider than those for k = 0. Hence, the projections are larger than those of the 
ground state tubes, which is clearly reflected by the wave functions. Again, it is possible to 
form a superposition that is localized along only one of the tubes (see Fig. 11), but now the 
tunneling time is 11 orders of magnitude shorter than that for k = 0. 
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The vortex tubes for k = 6, H e ff = 0.66219 and f3 = 0.5 still roughly "fit inside" the 
N = 2 islands of regular motion, see Fig. 1. But the quasienergy spectrum displayed in 
Fig. 4 indicates that the number of states accessible to the Mathieu approximation is larger 
than the area of the regular islands, divided by 27rft e //; the hierarchy of near- resonant states 
extends farther than one might expect. Further evidence for this statement is provided by the 
following two figures. Fig. 12 shows the probability density of a Floquet state with k = 10, 
whereas Fig. 13 shows the density of a superposition with its partner state. The tubes 
for k = 10 do not fit into the islands of regular motion; they would lie in the surrounding 
stochastic sea. Nevertheless, the wave functions are structured similarly to those in Figs. 10, 
11, as if there still were "quantizable" vortex tubes. Hence, although fully preserved vortex 
tubes do no longer exist outside the islands, it must be possible to use their remnants, such 
as Cantori, in a semiclassical quantization procedure. In this way, the quantum mechanical 
regularity associated with regular classical resonant motion is continued into a regime where 
a major part of the corresponding classical dynamics is chaotic. Expressed differently, the 
Mathieu quantum number k is not necessarily related to invariant vortex tubes; the number 
of Floquet states that can be labelled by k is larger than the number of vortex tubes that 
fit into the regular islands. Even though the transition from (mainly) regular to (mainly) 
stochastic classical motion appears to be quite sharp on the borders of the islands of stability, 
there is no correspondingly sharp transition in quantum mechanics. 



It had already been pointed out that the Mathieu equation (|23|) can be regarded as 
the Schrodinger equation for a fictitious particle in a cosine lattice. When the energy of 
such a particle is close to the top of the barrier, i.e., when afc(z/(j),g) ~ 2q, its probability 
density is strongly enhanced in the vicinity of the potential maxima, reflecting the fact that 
a corresponding classical particle would spend most of its time there. Just as the minima of 
the cosine potential correspond, within the pendulum approximation, to the stable, elliptic 
periodic orbits, the maxima correspond to the unstable, hyperbolic ones. According to 
this qualitative consideration, Floquet states with a Mathieu quantum number k such that 
Qk{v{j),Q) ~ 2g should be separatrix states, with a probability density that is strongly 
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enhanced along the unstable, hyperbolic periodic orbits. 

Quantitatively, the simple Mathieu approximation becomes less accurate for large quan- 
tum numbers k. For the present parameters, a characteristic value with k = 17 is clos- 
est to 2q, whereas it is actually a state with k — 16 whose density exhibits the clearest 
traces of the unstable classical orbit (see Fig. 14). Nervertheless, the approximation still 
yields the correct order of magnitude for the tunnel splitting: analytically, one obtains 
|(Ae)i 6 /^ e //| « 1.4 • 1CT 2 , compared to the numerically found value of 4.6 • 1CT 2 . Naturally, 
close to the top of the cosine wells the splitting is quite large (in the present example, the 
absolute value of (Ae)\Q is 21 orders of magnitude larger than the ground state splitting!), 
and the probability density of both members of a doublet shows significant differences. In 
spite of that, Fig. 15 demonstrates that it is still possible to form a superposition of a sep- 
aratrix state with its partner state such that, at least for a short time, the resulting density 
is strongly enhanced along only one of the two unstable 47r-periodic orbits. But complete 
destructive interference along the other one is no longer possible. The large tunnel splitting 
results in a quite short tunneling time for this doublet, r tunne i « 10.9 -27T. Thus, after merely 
11 cycles of the driving force the density will be concentrated along the other hyperbolic 
orbit. 

Quantum mechanical energy eigenstates with a density that peaks along unstable, hy- 
perbolic periodic orbits, first investigated by McDonald and Kaufmann and Heller ||45|| , 



have been dubbed "scars". In the context of periodically driven systems, such scars have 
attracted some interest since it was speculated that the experimentally observed anomalous 



stability of certain hydrogenic Rydberg states against microwave ionization [46] could be 



related to a scarred Floquet state [[47]]. In fact, the scar investigated by Jensen et al. |47| 



turns out to be simply the separatrix state of an N = 1 resonance. The present simple 
model shows, first of all, why an unstable classical orbit can "attract" quantum mechanical 
probability density, and it also shows under which conditions this phenomenon can occur 
in Floquet states. To get further insight into a possible connection between scarring and 
enhanced stability against ionization, it might be interesting to investigate also an iV = 2 res- 
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onance in highly excited, microwave-driven hydrogen atoms: are there, again, anomalously 
stable states? 



VII. IMPLICATIONS 

The preceding investigation of Floquet states and quasienergies linked to a primary N = 2 
resonance has demonstrated that the incompatibility of the 27r-periodic boundary conditions 
imposed on the quantum mechanical Floquet states, and the 2 • 2tt periodicity of classical 
invariant manifolds, necessarily leads to the tunnel effect. Pictorially speaking, there is 
quantum tunneling from one of the N = 2 islands seen in Fig. 1 through the stochastic sea 
to the other. Thus, one is faced with a special type of "dynamical tunneling", as studied 



by Davis and Heller [f|l| in the context of energy eigenstates. The generalization to the 
case of an arbitrary N is obvious: the Mathieu approximation leads to a cosine lattice with 
periodic boundary conditions after N wells. The wave functions associated with classical, 
N • 27r-periodic vortex tubes correspond to the single-well "atomic" states of that lattice, 
whereas the actual near-resonant Floquet states correspond to delocalized Bloch waves. 

Besides dynamical tunneling, another feature known from the investigation of time- 
independent systems carries over to the periodically time-dependent case. Jaffe, Shirts and 
Reinhardt |49|j50[| have argued that remnants of destroyed invariant manifolds, so-called 
"vague tori" existing in the apparently chaotic part of phase space, may introduce enough 
regularity to allow EBK-like quantization even in the absence of complete tori. Analogously, 
"vague vortex tubes" should be responsible for the surprisingly regular appearance of the 
Floquet state displayed in Fig. 12. The notion of a vague torus, although somewhat hard to 
formulate in precise mathematical terms, provides a bridge between regular and fully chaotic 
classical dynamics, and gives some insight how a quantum system behaves when its classical 
counterpart falls into this transition regime. The numerical example of the previous section 
also shows that there can be a closed link, expressed by the existence of a common quantum 
number k, between Floquet states associated with invariant manifolds surrounding stable, 
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elliptic periodic orbits and states associated with unstable, hyperbolic periodic orbits. In the 
general case, one might wish to apply torus quantization to the classically regular, elliptic 



motion and a periodic orbit cycle expansion to the chaotic, homoclinic motion ||13|| . But to 
which extent, if at all, is it possible to "interpolate" between both quantization schemes? 
One-dimensional, periodically time-dependent systems like the one studied here, with their 
merely three-dimensional phase space {(p, x, r)}, might provide the simplest access to this 
important problem. 

Returning now to the specific properties of the triangular well model, it has been sug- 
gested that there exists a transition from a pure point quasienergy spectrum to a continuous 



one at a nonzero value of the amplitude (3 [p8| , p9| . A point spectrum would be connected 



with square-integrable Floquet states that have a "localized" distribution of coefficients 
when expanded in the basis of unperturbed triangular well eigenstates; a continuous spec- 
trum would be connected with "extended" states. The classical dynamics indeed shows a 
qualitative change when f3 is enlarged: for small amplitudes, all trajectories are either tied 
to invariant vortex tubes or move in the stochastic layers surrounding the individual reso- 
nance zones, but remain bounded. Within a comparatively small interval of (3, all resonances 
overlap (see fllTD), so that the individual stochastic layers merge into a connected stochastic 
sea. We then have a dichotomy: the trajectories are either confined to the remnants of the 
regular resonant islands or they move in the stochastic sea (cf. Fig. 2), where they can gain 
an unlimited amount of energy. It would be tempting to speculate that the latter type of 
behavior is linked to "extended" Floquet states in the quantum system. 

The numerical evidence remains inconclusive. Fig. 16 shows the near-resonant 
quasienergy spectrum (N = 2) for h e ff = 0.66219, for amplitudes (3 between and 1.0. 
With increasing (3, more and more resonant vortex tubes are destroyed (cf. Figs. 1 and 2), 
and the number of states that can be described by the Mathieu approximation is diminished. 
Nevertheless, there are clear signs of the resonance-induced regularity even for (3 = 1.0. 
Fig. 17 shows the probability density of a member of the resonant ground state doublet for 
this amplitude, whereas Fig. 18 shows the density of a "one-tube" superposition with its 
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partner state. Although these wave functions appear considerably more "ragged" than they 
did for (3 = 0.5 (cf. Figs. 8 and 9), there does not seem to be a sharp transition, but merely 
a gradual change of behavior. However, the quantum mechanics of the periodically driven 
triangular well is more subtle than a numerical investigation might suggest. 

The previous analytical deliberations were concerned with single resonances. But for 
a rigorous study of the nature of the quasienergy spectrum the single-resonance approxi- 
mation is insufficient. To get a glimpse of the problem, consider the approximate Mathieu 



spectrum ( 52|) for the iV-th resonance, together with that for the (N + l)-th, projected into 
one Brillouin zone. At certain values of j3, approximate eigenvalues belonging to the iV-th 
resonance will cross others of the (N + l)-th. However, the full Hamiltonian of the driven 
triangular well has no symmetry which could allow these crossings: the terms neglected in 
the derivation of the two near-resonant spectra will lead to a small, but nonvanishing cou- 
pling of the Floquet states belonging to the two different resonances. Hence, the eigenvalues 
will repel each other [[Hj], so that the crossing of the approximate eigenvalues become anti- 



crossings, most of them much too tiny to be detectable by numerical means. Accordingly, 
the Floquet states become linear combinations of the approximate single-resonance states. 
In other words, there is quantum tunneling not only between invariant manifolds belonging 
to the same N, but also between those belonging to different N. Now there are not only two 
primary resonances, but infinitely many. According to fll6|) , the width Aijv °f the classical 
resonance zones grows quadratically with N, so that even for arbitrarily small f3 there will 
always be an infinite number of primary resonances that support an arbitrarily large number 
of Floquet states, no matter how large (or small) the effective Planck constant h e ff- If one 
imagines all the infinitely many near-resonant quasienergy spectra superimposed in a single 
Brillouin zone for, say, < (3 < 1, with all crossings replaced by anticrossings, one gets an 
extremely complicated net, with arbitrarily narrow anticrossings appearing everywhere. It 
seems unavoidable to conclude that the true quasienergies of the driven triangular well can 
not be differentiable with respect to the amplitude f3. 

A question that could shed some light on the nature of the quasienergy spectrum can 
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therefore be phrased as follows: if one starts from the approximate near-resonant Floquet 
states, and then takes into account quantum tunneling between different resonances, are the 
resulting Floquet states "localized" over a finite number of resonances, or "extended"? It 
appears quite possible that, because of tunneling, there is no spectral transition except for 
(3 = 0, and all Floquet states might be extended for every (3 > 0. Or could localized and 
extended states coexist, such that the number of extended states gradually increases when 
(3 is enlarged? 

For a physicist, the approximate analytical and numerical results provide a description 
of the dynamics over time scales that are short compared to inter-resonance tunneling times, 
which may be sufficient for most practical purposes. From the point of view of a rigorous 
mathematician, the present study of the periodically driven triangular well certainly poses 
additional questions. It suggests, however, a more direct starting point for the investigation 
of periodically time dependent quantum systems than the one adopted in p0|-p5|. The 
model of the kicked rotor |TJ| , which was motivated by a single- resonance approximation, has 
meanwhile been understood in some detail [r5, |T5| i P0|,pi| . The next generation of problems 



seems to arise from the interaction of different resonances; a detailed analysis of the role 
of quantum tunneling is required. The triangular well model P?^p5| appears to be ideally 
suited for that purpose. 
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FIGURES 

FIG. 1. Poincare surface of section for the driven particle in a triangular well potential (0), 
taken at r = 2>tt/2 mod 2ir. The island in the lower left corner originates from the N = 1 resonance, 
the two large islands from the resonance with N = 2. The driving amplitude is (3 = 0.5, lower 

(c) 

than p\ = 0.537 which is the value required by the Chirikov criterion for the overlap of these two 
resonances. The two boxes in the upper right corner enclose areas 2irh e ff, with h e fj = 1.33325 
and 0.66219, respectively. 

FIG. 2. Poincare section for (3 = 1.0. The driving amplitude now is significantly higher than 

(c) 

Poo = 0.617, above which, according to the Chirikov criterion, all resonances overlap. Nevertheless, 
there are still remnants of the regular islands for N = 1 and N = 2. The two boxes in the upper 
right corner have the same meaning as in Fig. 1. 

FIG. 3. Part of the numerically computed quasienergy spectrum for Ti e ff = 1.33325, restricted 
to one Brillouin zone. Energies of unperturbed triangular well eigenstates close to the "resonant" 
state r = 20 differ by approximately H e ff/2. The behavior of the quasienergies of near-resonant 
states is well described by the Mathieu approximation (|32|); integers in the right margin are the 
quantum numbers k for one of the two resonant ladders. 

FIG. 4. Part of the numerically computed quasienergy spectrum for h e ff = 0.66219. For this 
value of h e ff, energies of unperturbed triangular well eigenstates close to the "resonant" state 
r = 40 differ by approximately fi e ff/2. 

FIG. 5. Approximate quasienergy spectrum according to the Mathieu approximation (|52|), for 
the same parameters as used in Fig. 4, and k = 0, . . . , 9. 

FIG. 6. Negative dipole matrix elements (n|x|n + 2) in the vicinity of n = 40. Within the 
Mathieu approximation leading to Fig. 5, it is assumed that all these elements are equal to unity. 
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FIG. 7. Contour plot of the probability density of a member of the ground state doublet (k = 0), 
for (3 = 0.5 and H e ff = 1.33325. The density is concentrated along the two classically equivalent 
"innermost quantized vortex tubes" surrounding the two central elliptic periodic orbits of the 
N = 2 resonance, projected to the (x, r) plane. Strong interference occurs where both projections 
intersect. In this and the following figures, the displayed interval of space ranges from x m i n = to 
Xmax = 25, the interval of time corresponds to 4 • 2ir, i.e., to 4 cycles of the driving force. 

FIG. 8. Probability density of a member of the ground state doublet (k = 0), for (3 = 0.5 and 
H e ff = 0.66219. The wave function is now tied closer to the associated classical object than it is 
in the case of the larger h e ff, see Fig. 7. 

FIG. 9. Probability density of a solution to the time-dependent Schrodinger equation that 
consists of a superposition of both members of the ground state doublet. The parameters are again 
(3 = 0.5, H e ff = 0.66219. The interval of time corresponds to 4 cycles of the driving force. The 
wave function appears to follow only one of the two quantized vortex tubes, and thus exhibits 
subharmonic motion. However, it is coupled to the other tube by the tunnel effect. The tunneling 
time is r tunneX « 4 • 10 22 • 2vr. 

FIG. 10. Probability density of a member of the doublet with k = 6. The associated vortex 
tubes still roughly fit inside the regular elliptic islands, cf. Fig. 1 ([3 = 0.5, h e jf = 0.66219). 

FIG. 11. Solution to the time-dependent Schrodinger equation that consists of a "one-tube" 
superposition of the k = 6 doublet ((3 = 0.5, h e jf = 0.66219). 

FIG. 12. Probability density of a member of the doublet with k = 10. Although the associated 
vortex tubes would not fit inside the regular elliptic islands seen in Fig. 1, this eigenfunction is 
structured analogously to that displayed in Fig. 10 ((3 = 0.5, fi e ff = 0.66219). 

FIG. 13. Solution to the time-dependent Schrodinger equation that consists of a superposition 
of both members of the k = 10 doublet {(3 = 0.5, Ti eff = 0.66219). 
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FIG. 14. Probability density of a member of the k = 16 doublet. This is a separatrix state: its 
density is concentrated along the projections of the two hyperbolic periodic orbits originating from 
the N = 2 resonance {(3 = 0.5, h e ff = 0.66219). The arrows on the bottom indicate the moments 
where the stable, elliptic periodic orbits hit the wall, cf. Fig. 7. 

FIG. 15. Solution to the time-dependent Schrodinger equation that consists of both members 
of the k = 16 doublet. This wave function shows signs of "scarring", i.e., its density appears 
strongly enhanced along the projection of a single hyperbolic periodic orbit of the N = 2 resonance. 
However, the tunneling time for this doublet is merely 10.9 -2tt, so that the density will be enhanced 
along the other hyperbolic periodic orbit after merely 11 cycles (/? = 0.5, fi e ff = 0.66219). 

FIG. 16. Part of the quasienergy spectrum for h e ff = 0.66219, and amplitudes (3 ranging from 
to 1.0. Even for the highest amplitude, there are clear signs of the resonance-induced regularity. 

FIG. 17. Probability density of a Floquet state with k = for h eff = 0.66219 and = 1.0, 
cf. Fig. 2. 

FIG. 18. Superposition of both members of the N = 2 ground state doublet for h e ff = 0.66219 
and = 1.0. 
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